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1. Introduction 

Heavy quarkonia can play an important role in studying hot and dense strongly interacting 
matter. Because of the heavy quark mass quarkonia binding can be understood in terms of the 
static potential. Also becasue of their relatively small size they do not experience the medium till 
relatively high temperatures. General considerations suggest that quarkonia will melt at tempera- 
tures right above the deconfinement temperature as a result of modification of inter-quark forces 
(color screening). Thus it was conjectured by Matsui and Satz [p that melting of different quarko- 
nia states due to color screening can signal Quark Gluon Plasma formation in heavy ion collisions. 
Classic studies of quarkonia dissolution [Q, ^, ^ ^, ^, ^ rely heavily on potential models. However 
it is very unclear if such models are valid at finite temperature As an alternative method to 
approach this problem we may use lattice to determine meson spectral functions. For charmonium 
such studies appeared only recently and suggested, contrary to potential models, that and T]c 
survive at temperatures as high as 1.67^ [Q, [l^ 11]. It has been also found that Xc melts at tem- 



perature of about 1.1 Tc [|11|]. In principle, meson spectral functions offer an attractive possibility 
to study the properties of radially excited quarkonium states. In particular it is straightforward to 
determine quai^konium decay width to dilepton. In recent studies of charmonium spectral functions 



several peaks were identified beyond and ground state M, hQ, 11]. It is not clear, however, if any 



of those peak correspond to radial excitation. Moreover it has been argued that all peaks in the 



spectral functions except the first one are lattice artifacts [|11[]. 

In this contribution we focus on our updated results for the charmonium; bottonium results 
and detailed study of the Maximum Enthropy Method reliability will be presented in our upcoming 
paper. Preliminai'y results for bottomonium can be also found in Ref. 



2. Meson correlators and spectral functions 

In our lattice investigation we calculate correlators of point meson operators of the form 

JH{t,x)=q{t,x)THq{t,x), (2.1) 

where Tu = l,75,7^,757p,7/i7v fixes the quantum number of the channel to scalar-, pseudo-scalar, 
vector, axial-vector and tensor channels correspondingly. The relation of these quantum number 
channels to different meson states is given in Tab. [1} 



r 








cc{n = 1) 


cc{n = 2) 


75 




0-+ 




Tic 


1 


7- 




1 




J/v 


V 


77' 




1+- 




he 




1 


'Po 


0++ 




XcO 




757. 


'Pi 


1++ 




Xcl 








2++ 




Xc2 





Table 1: Charmonia states in different channels 
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Most dynamic properties of a finite temperature system are incorporated in the spectral func- 
tion. The spectral function Oh{pq,p) for a given mesonic channel // in a system at temperature T 
can be defined through the Fourier transform of the real time two-point functions and D< or, 



equivalently, as the imaginary part of the Fourier transformed retarded correlation function []130, 

Oh{po,p) = ^iD^{Po,p)-D^{Po,p)) = (2.2) 



-ImDfi{po,p) 



K 

4 



D^h' '(PoJ) = I j^.e'P'^D^^'^'ixo^x) (2.3) 



D^(xo,x) = {Jh{xq,x),Jh{0,0)) 

D^{xo,x) = {JHiO,0), Jh {xq^x)), xq>Q (2.4) 



The Euclidean time correlator calculated on the lattice 



Gh{^,p) = J d\e'P\T,JH{T,x)JH{0M (2.5) 

is an analytic continuation of the real time correlator Ge{'^,p) = {—iz,p). 

Using this equation and the Kubo-Martin-Schwinger (KMS) condition [|l3|] for the correlators 

D^{xo,x)=D<{xo + i/T,x), (2.6) 



one can relate the Euclidean propagator Gh{'T:,p) to the spectral function, Eq. (2.2), through the 
integral representation 



poo 

G{t,p) = / d(oo{(0,p)K{(o,T) 
Jo 



cosh(oj(T-l/2r)) 

K{(0,t) = . . ■ (2-7) 

smh(ft)/2r) 

To reconstruct the spectral function from the lattice correlator G{z,T) this integral represen- 
tation should be inverted. Since the number of data points is less than the number of degrees of 
freedom (which is ^(100) for reasonable discretization of the integral ) spectral functions can be 



reconstructed only using the Maximum Entropy Method (MEM) [14|. In order to have sufficient 
number of data points either very fine isotropic lattices [|ri|] or anisotropic lattices |^ |T^ should 
be used. Another difficulty arises due to the large quark mass; discretization eiTors 0'{amc^b) are 
present in the heavy quai^k system. To remove these discretization eiTors we use the Fermilab 
approach [[l^, [l6|]. 

We performed calculation of chaimonia correlators in quenched QCD using anisotropic lat- 
tices and Wilson gauge action. For charmonium we use following gauge coupling values jS = 
5.6, 5.7, 5.9, 6.1 with anisotropy = av/a, = 2 and /3 =6.1, 6.5 with ^ = 4. When the lattice 
spacing is set using the Sommer scale ro = 0.5 fm the above gauge coupling con^espond to lat- 
tice spacing a^^ = 1.56, 1.91, 2.91, 4.11, 8.18, 14.12 GeV. Further details of the lattice action 
together with the parameters can be found in Ref. [[l^, p^ ]. 

To be able to reconstruct the spectral functions high statistical accuracy for charmonia coitc- 
lators is required. This makes calculations computationally intensive, so we performed them on a 
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prototypes of RBRC QCDOC machine. It is a dedicated Lattice QCD machine developed by physi- 
cists from Columbia University, BNL, RIKEN and UKQCD. Three such machines, each reaching 
about lOTFlops peak performance, are currently installed at BNL and EPCC. Our prototypes were 
single-motherboard machine at about 50 GFlops peak. Such resources are still not adequate for the 
full QCD simulations, so we use quenched approximation, which is equivalent to neglecting quark 
loops. Typical statistics gathered was 500 to 1000 measurements, separated by 400 updates. 



3. Bayesian analysis of meson correlators 



In Bayesian analysis of the correlator one looks for a spectral function which maximizes the 
conditional probability P[o\DH] of having the spectral function a given the data D and some prior 

). Different Bayesian methods differ in the choice of the 



knowledge H (for a reviews see [ ]14[ , |]J 
prior knowledge. One version of this analysis which is extensively used in the literature is the 
Maximum Entropy Method (MEM) [12, |l^]. In this method the basic prior knowledge 

is the positivity of the spectral function and the prior knowledge is given by the Shannon - Janes 
entropy 



d(0 



a(w)-m(w)-a(w)ln(^-^) 

m ft) 



(3.1) 



The real function m{co) is called the default model and parametrizes all additional prior knowledge. 



about the spectral functions, such as the asymptotic behavior at high energy [14]. For this case the 
conditional probability 



P[a\DH] = txp{ 



(3.2) 



with being the standard likelihood function and a a real parameter. To maximize P[o\DH] 



usually the Bryan algorithm [ |20| ] is used. We use a new approach (to be described in []17[]) which 
through exact mathematical transforamtions reduces the problem to the task of minimization of a 
function with positive definite second derivative. This is a minimization problem in number-of-data 
dimensions, which is not any more difficult task than applying a method. The Bryan algorithm 
takes different approach of singular value decomposition in order to find a relevant subspace, where 
minimization is performed. 



4. Spectral functions at zero temperature 

The spectral function for pseudo-scalar charmonium spectral functions is shown in Fig. 1. 
The first peak in the spectral function corresponds to r\c{^S) state. The position of the peak and 
the corresponding amplitude (i.e. the area unders the peak) are in good agreement With the results 
of simple exponential fit. The second peak in the spectrak function is most likely the combination 
of the excited states as its position and amplitude is higher than what one would expect for 2S 
state. The spectral function becomes sensitive to the effects of the lattice spacings for (O > 5GeV. 
In this (O region the spectral functions becomes sensitive to the choice of the default model. This is 
because only a very few data points in the con^elator carry information about the spectral function 
in that region. 
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Figure 1: Charmonium spectral function in the psuedoscalar channel (left) and the scalar channel (right) at 
different lattice spacings and zero temperature. 



Also shown in Fig.l is the spectral function in the scalar channel. The 1st peak corresponds to 
Xco{^P) state. The correlator is more noisy in the scalar channel than in the pseudo-scalar one. As 
the results the Xc{){^P) peak is less pronounced and has larger statistical errors. The peak position 
and the area unders the peak is consistent with the simple exponential fit. As in the psudo-scalar 
case individula excited states are not resolved and the spectral function depends on the lattice 
spacing and default model for CO > 5GeV. 

Similar results have been found for the vector and axial-vector channels which correspond to 
J /\lf and Xc\ states respectively. 

We would like to know what happens to different charmonia states at temperatures above the 
deconfinement temperature T^. With increasing temperature it becomes more and more difficult to 
reconstruct the spectral functions as both the number of available data points as well as the physical 
extent of the time direction (which is \/T) decreases. Therefore it is useful to study the temper- 
ature dependence of charmonia coiTclators first. From Eq. ( ^77| ) it is clear that the temperature 
dependence of charmonia coiTclators come from two sources: the temperature dependence of the 
spectral function and temperature dependence of the integration kernel K{t,0),T). To separate 
out the trivial temperature dependence due to the integration kernel, following Ref [ jlTl ] at each 
temperature we calculate the so-called reconstructed correlator 

poo 

Grecon{T,T)= dcoo {co,T = 0)K{r , (0,T) (4.1) 
JO 

Now if we assume that there is no temperature dependence in the spectral function - then the ratio of 
the original and the reconstructed correlator should be close to one, G{z,T)/ Grecon {''^^T) This 
way we can identify the cases when spectral function itself changes dramatically with temperature. 
This gives reliable information about the fate of charmonia states above deconfinement. 

In Fig. ^ we show this ratio at j8 = 6.5 and 6.1 for pseudoscalar and scalar channels of char- 
monium correspondingly. 

>From the figures one can see that the pseudo-scalar correlators shows only very small changes 
till l.STc indicating that the r\c states survives till this temperature with little modification of its 
properties. On the other hand the scalar correlator shows large changes already at 1 . 167^ suugesting 
strong modification or dissolution of the XcO state at this temperature. 
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Figure 2: The ratio G{T,T)/Grecim{'^,T) of charmonium for pseudoscalar channel at at a, ^ = 14.11GeV 
(left) and scalar channel at at a^^ ~ 8.18GeV (right) at different temperatures. 
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Figure 3: Charmonium spectral function in the psuedoscalar channel at a^^ = 14. 1 IGeV (left) and the scalar 
channel (right) at a,^^ ~ 8.18GeV at zero and above deconfinement temperatures. For finite temperature 
scalar chanel two different default models are shown. 



More detailed information on different charmonia states at finite temperature can be obtained 
by calculating spectral functions using MEM. The result of these calculation is show in Figs. ^. 
Because at high temperature the temporal extent and the number of data points where the coiTclators 
are calculated become smaller the spectral functuon reconstructed using MEM are less reliable. To 
study the temperature modifications of the spectral function we compare the finite temperature 
spectral functions against the zero spectral functions obatined from the correlator using the same 
time interval and number of data points available at finite temperature. We see that spectral function 
in the pseudo-scalar channel show no temperature dependence within the statistical errors shown 
in the figure in accord with the analysis of the correlation functions. Also the spectral functions 
show very little dependence on the default model. 

The scalar spectral function on the other hand shows large changes at 1.167^ which is consis- 
tent with coiTclator-based analysis. Also default model dependence of the scalar correlator is large 
above the deconfinement transition (c.f. Fig. 3, right). 
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